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2.0  TECHNICAL  DISCUSSION 


2 . 1  BACKGROUND 

2.1.1  Background  Summary 

The  basic  system  operation  is  explained  in  Ref.  6.  A  brief 
summary  of  the  system  operation  is  as  follows.  Consider  a  sample 
consisting  of  a  single-layer  thin  metal  coating  deposited  on  a 
transparent  dielectric  substrate.  The  objective  is  to  measure 
the  thermal  transport  properties  of  the  coating/substrate  system. 
The  method  of  measurement  is  a  new  technique  called 
d 1 f f us i on - wa ve  interferometry.  The  measurement  involves  two 
steps:  sinusoidal  heating  of  the  metal  coating  surface,  and 

measurement  of  the  sinusoidal  temperature  response  of  the  metal 
coating  surface. 

2.1.2  Theoretical  Background 

Under  one-d i mens i ona  1  operating  conditions,  the  sinusoidal 
heating  generates  a  critically-damped  steady  state  temperature 
field  known  as  a  thermal  diffusion  wave: 

T(x,t)  =  Tq  e  cos(kx  -  ojt  +  $ ) 
where  T  is  temperature,  k  is  wavenumber,  oj  is  angular  frequency. 

X  is  distance  and  t  is  time  (ref.  6).  Because  the 
thermal-diffusion  wave  is  partially  reflected  at  the 
coating-substrate  interface,  the  sinusoidal  temperature  response 


at  the  metal  coating  surface  is  out  of  phase  with  the  sinusoidal 
heating.  Changing  the  modulation  frequency  of  the  heating  causes 
the  wave  number  k  of  the  thermal  diffusion  wave  to  change. 
Measurement  of  the  phase  shift  between  the  heating  and  the 
temperature  response  as  a  function  of  the  modulation  frequency 
yields  an  interference  oscillation  in  the  phase  that  allows  the 
thermal  transport  properties  of  the  coating/ substrate  system  to 
be  determined. 

2.1.3  Experimental  Background 

The  sinusoidal  heating  is  produced  using  an  argon-ion  pump  laser 
whose  intensity  is  modulated  by  an  electrooptic  modulator.  The 
temperature  response  is  measured  by  reflecting  an  unmodulated 
HeNe  laser  beam  off  the  sample,  and  measuring  the  intensity 
modulation  induced  by  sinusoidal  modulation  of  the  metal  coating 
reflectivity  (caused  by  the  temperature  dependence  of  the  metal 
coating  reflectivity).  The  reflected  HeNe  beam  is  monitored 
using  a  photoconduc t  i  ve  detector,  and  a  lock-in  analyzer  is  used 
to  measure  the  resulting  weak  synchronous  signals. 

2 . 2  SYSTEM  UPGRADE 

2.2.1  Pump  Laser 

The  original  system  included  a  100  mW  argon-ion  pump  laser.  This 
laser  was  adequate  for  the  initial  studies.  However,  the  low 
power  required  that  a  fairly  small  spot  diameter  he  used  in  order 
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to  obtain  sufficient  heat  l hr  to  yield  an  acceptable  signal  to 
noise  ratio.  For  simplicity  of  data  analysis,  it  is  desirable 
that  radial  heat  flow  be  negligible  at  all  modulation  frequencies 
employed.  This  one-dimensional  operating  requirement  is 
equivalent  to  requiring  that  the  thermal  diffusion  length  be  much 
less  than  the  pump  spot  radius  at  all  modulation  frequencies 
employed.  Therefore,  the  lowest  operating  frequency  within  the 
one-dimensional  limit  is  limited  by  the  spot  diameter  used.  In 
order  to  increase  the  system  dynamic  range  over  which  the 
modulation  frequency  may  be  varied  within  one-dimensional 
operation,  the  100  mW  pump  laser  was  replaced  by  6W  argon-ion 
laser.  As  will  be  shown  later,  attaining  the  greatest  possible 
dynamic  range  is  the  key  to  assuring  uniqueness  of  parameter  fit 
when  analyzing  the  data. 

2.2.2  Electro optic  Modulator 

An  electrooptic  modulator  Is  used  to  modulate  the  intensity  of 
the  pump  laser  beam.  Although  the  original  electrooptic 
modulator  was  specified  to  operate  at  laser  powers  up  to  10W.  it 
was  found  that  modulated  beam  quality  was  unacceptable  above 
100  mW  of  pump  power,  apparently  due  to  thermal  lensing  in  the 
index  matching  fluid.  A  custom-designed  electrooptic  modulator 
with  reduced  optical  path  length  within  the  index  matching  fluid 
was  obtained  under  warranty  from  the  modulator  manufacturer. 
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2.2.3  Signal  Detector 

A  photoconductlve  detector  is  used  to  detect  the  weak  intensity 
modulation  induced  on  the  probe  HeNe  beam.  The  detector  housing 
originally  in  use  failed  to  provide  adequate  shielding  at 
frequencies  above  10  kHz.  The  temporary  solution  to  this  problem 
was  to  measure  the  synchronous  background,  and  subtract  it  from 
the  signal  plus  background.  This  subtraction  resulted  in  a 
reduction  of  the  signal  to  noise  ratio  by  a  factor  of  about  2. 

An  improved  detector  housing  was  installed,  reducing  the 
synchronous  background  to  a  negligible  level. 

2.2.4  Rear-surface  Probe 

The  system  configuration  at  the  start  of  the  contract  allowed  the 
temperature  modulation  of  the  film/air  interface  to  be  monitored 
using  a  HeNe  laser  relying  on  the  temperature  modulation  of  the 
metal  coating  reflectivity.  An  additional  probe  layout  has  been 
installed  to  allow  temperature  modulation  measurements  at  the 
f i 1  m/subs t ra te  interface  for  metal  coatings  deposited  on  a 
transparent  substrate  (See  kef.  5  for  a  similar  technique).  In 
this  case,  a  HeNe  laser  beam  passes  through  the  substrate  to 
strike  the  metal  coating  at  a  point  centered  within  the  region 
heated  by  the  pump  laser  beam.  After  being  reflected  by  the 
metal  coating,  the  rear  II »•%’*•  probe  beam  exits  through  the 
substrate,  and  after  passing  through  a  relay  lens  enters  a  • 
photoconductlve  detector  identical  to  that  used  to  monitor  the 
front-surface  probe  beam.  The  relay  lens  is  needed  to  form  an 


Imaging  system.  Without  the  relay  lens,  probe  beam  deflection 
introduced  by  a  pump-induced  thermal  lens  in  the  substrate  would 
otherwise  give  rise  to  a  strong  synchronous  signal  due  to 
detector  reponse  inhomogeneity.  Data  from  the  rear-surface  probe 
provides  at  least  one  independent  condition  when  fitting 
adjustable  parameters  during  data  anlysis.  This  addtional 
independent  data  correspondingly  reduces  the  necessary  system 
dynamic  range  required  to  insure  uniqueness  of  parameter  fit. 

2 • 3  SYSTEMATIC  ERROR 

2.3.1  Details 

A  major  systematic  error  has  been  discovered  and  eliminated.  The 
error  entered  during  the  calibration  of  system  electronic  phase 
shift.  As  discussed  in  section  2.1.  the  phase  shift  between  the 
sinusoidal  heating  of  the  metal  coating  surface  and  the 
sinusoidal  temperature  response  of  the  metal  coating  surface 
constitutes  the  desired  data.  The  lock-in  analyzer  provides  a 
phase  shift  value,  but  this  value  Includes  electronic  phase  shift 
introduced  In  the  detector  and  preamplifier.  In  order  to 
calibrate  this  electronic  phase  :  ,ift.  the  pump  beam  Is  directed 
Into  the  probe  detector.  A  nonzero  phase  shift  is  assumed  to  be 
electronic  in  origin.  However,  an  error  occurred  due  to  the 
presence  of  a  red  colored  glass  filter.  This  filter  passes  HeNe 
laser  light,  but  absorbs  argon  laser  light.  Tin*  filter  is  used 
to  eliminate  scattered  pump  light  during  probe  beam  measurements. 


However,  during  calibration  of  the  electronic  phase  shift,  the 
pump  bean  was  originally  directed  through  the  glass  filter  into 
the  probe  detector.  Precautions  were  taken  to  reduce  the 
intensity  of  the  pump  light  sufficiently  to  reduce  heating  in  the 
glass  filter  to  a  (hopefully)  negligible  level.  The  criterion 
used  was  that  the  measured  phase  be  independent  of  the  pump  power 
over  a  two  decade  range  of  pump  power.  Unfortunately,  this 
criterion  proved  to  be  inadequate,  and  a  large  difference  in 
phase  was  found  to  occur  between  calibration  readings  taken  with 
and  without  the  glass  filter. 

2.3.2  I  up  1 Ications 

The  discovery  of  a  major  systematic  error  after  two  years  of 
operation  serves  to  Illustrate  the  need  for  a  wide  dynamic  range 
of  modulation  frequency  and  consequently  of  thermal 
d  i  f f us  ion -wave  length.  Had  this  error  not  been  discovered  by 
accident,  a  problem  would  have  been  ultimately  encountered  in 
fitting  the  data  over  the  new  1  y - i nc r e a s ed  dynamic  range. 

However,  if  would  not  have  been  clear  whether  the  problem  was  due 
to  the  presence  of  a  systematic  error,  or  to  the  inadequacy  of 
the  sample  model,  say  due  to  supposed  coating  inhomogeneity. 

Many  such  sizeable  systematic  errors  have  already  been  discovered 
and  eliminated.  More  such  errors  may  still  be  present. 
Measurements  will  continue  on  the  simplest  possible  samples  (i.e. 
single  metal  layers,  or  solid  metal  mirrors)  until  there  is  good 
cause  to  believe  that  the  measurements  are  free  of  uon-negl igible 
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systematic  errors. 


2 . 4  RESULTS 

2.4.1  Description  of  Figures 

Data  have  been  measured  on  a  1000  %  Ni  coating  deposited  on  an 
S102  substrate  by  e-beam  deposition.  The  data  are  shown  in  Figs. 
1  and  2,  and  representative  theoretical  model  behavior  is  shown 
in  Figs.  3  and  4.  for  comparison.  (The  theoretical  plots  are  not 
the  result  of  fitting  the  experimental  data.)  Figs.  1  and  3  show 
the  phase  of  the  sinusoidal  temperature  response  of  the  coating 
at  the  front  surface  (air/coating  interface)  relative  to  the 
sinusoidal  heating  by  a  pump  beam  incident  on  the  front  surface. 
Figures  2  and  4  show  the  phase  of  the  sinusoidal  temperature 
response  of  the  coating  at  the  rear  surface  (coating/substrate 
interface)  relative  to  the  sinusoidal  heating  by  a  pump  beam 
incident  on  the  front  surface.  The  behavior  is  plotted  as  a 
function  of  the  modulation  frequency  in  all  figures  (In  Figs.  3 
and  4,  x  is  proportional  to  the  square  root  of  the  modulation 
frequency).  Three  curves  are  shown  in  Figs.  3  and  4.  They 
correspond  to  different  relative  thermal  offusivities  of  the 
coating  and  substrate.  The  g  =  10  curves  correspond  to  a 
thermally  insulating  coating  on  a  thermally  conductive  substrate. 
The  g  *  0.1  curves  correspond  to  a  thermally  conductive  coating 
on  a  thermally  insulating  substrate.  The  g=l  curves  correspond 
to  thermally  matched  coating  and  substrate.  The  plotted  case  of 
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Fig.  2  Rear-Surface  Phase  vs  Modulation  Frequency 
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log  x [arbitrary  unitsj 


h  -  100  corresponds  to  a  relatively  large  thermal  contact 
resistance. 

2.4.2  Detailed  Theoretical  Behavior 

In  Figs.  3  and  4  the  oscillations  are  the  result  of  interference 
between  the  direct  heating,  and  the  diffusion  waves  that  are 
partially  reflected  at  interfaces.  The  interference  oscillations 
in  Fig.  3  are  confined  to  the  range  -45°  <  <  45°  where  ^  is 

phase.  The  oscillations  decrease  exponentially  in  amplitude  for 
arguments  greater  than  those  plotted  because  of  the  critically 
damped  nature  of  diffusion  waves.  The  lowest  frequency  zero 
crossing  shown  is  present  only  for  non-negl lglble  values  of  the 
thermal  contact  resistance  between  the  coating  and  the  substrate. 
The  phase  In  Fig.  4  decreases  mono  ton i ca 1 1 y  for  arguments  greater 
than  those  plotted  because  of  the  finite  transit  time  for  thermal 
waves  to  propagate  through  the  coating.  (The  diffusion  waves  are 
generated  at  the  front  surface  but  detected  at  the  rear  surface 
In  this  case . ) 

2.4.3  Discussion  o  f  Experiment  a  1  l)a_t  a_ 

The  zero  crossings  in  both  figures  Imply  t. he  presence  of  a 
non-negl lglble  contact  resistance.  The  two  data  sets  plotted  in 
Fig.  1  were  measured  about  one  week  apart  at  approximately  the 
same  point  on  the  sample.  The  source  of  the  nonreproduc i bi 1  i  t y 
of  the  low-frequency  data  is  not  known;  additional  studies  are  in 
progress  to  provide  an  explanation.  Although  a  three-parameter 
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model  has  been  derived  lo  fit  the  data,  a  three-parameter  fitting 
routine  has  not  been  implemented  yet.  This  joining  of  theory  and 
experiment  is  the  next  major  milestone  to  be  met.  At  the  present 
stage,  the  experimental  data  appear  to  be  consistent  with  the 
derived  model  behavior  (within  the  general  trends  that  can  be 
identified  without  knowing  the  values  of  the  fitting  parameters 
corresponding  to  the  experimental  data). 


2 . 5  THEORETICAL  MODELLING 

2.5.1  Background 

One  of  the  goals  of  this  research  is  to  characterize  the  thermal 
transport  properties  of  dielectric  coatings.  In  -  order  that  the 
thermal  diffusion  wave  be  generated  in  a  layer  thin  compared  to 
the  coating  thickness,  it  is  necessary  to  deposit  a  thin  metal 
layer  over  a  dielctric  coating  under  study.  The  presence  of  a 
second  coating  layer  introduces  additional  interference 
structure  Into  the  data,  which  must  be  modelled.  An  analytic 
boundary  value  solution  to  the  thermal  diffusion  equation  has 
been  obtained  under  the  following  assumptions:  two  coating 
layers,  sinusoidal  heating,  one -d i mens i ona 1  heat  flow  normal  to 
the  surface,  piecewise  constant  approximation,  neglect  of  the 
effect  of  nonzero  optical  depth  of  the  pump  laser  beam,  infinite 
substrate  thickness  and  thermally  thin  interfaces  of 
nonnegl igible  thermal  contact  resistance.  The  approach  is  a 
generalization  of  that  described  in  Ref.  6. 


12 


2.5.2  Ana  1 y a  1  s 


2 . 5 . 2 . 1  Boundary  Value  Equation 

The  temperature  and  heat  flux  at  the  surface  may  be  related  to 
the  temperature  and  heat  flux  at  the  substrate  by  the  following 
equation  ( Ref  .  6.8). 


where  T  -  temperature 
f  *  heat  flux 
=  -K  V  T 

where  Kj  *  thermal  conductivity  of  layer  j 

h  ■  SV'*" 

where  i  » 

kj  =  thermal  diffusion-wave  number 

-  ("f/Kj)172 

where  f  »  modulation  frequency 

ICj  *  thermal  diffusivlty  of  layer  j 
and  where  Cj  -  cosh  kj  ^ (  1  -  i  ) 

Sj  «  sinh  kj  j I  I  ►  1  ) 

where  *  thickness  of  layer  j 
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and  where  R^j  =  thermal  contact  resistance  of  unit  area 

between  regions  1  and  j. 

Subscripts  1,  2,  and  3  refer  to  the  outer  coating  layer  (metal), 
the  Intermediate  coating  layer,  and  the  substrate  respectively. 
This  model  can  be  used  to  describe  two  cases:  either  that  of  a 
dielectric  layer  with  a  metal  overcoat,  and  thermally  thin 
interfaces  between  the  layers,  or  that  of  a  metal  coating  with  a 
thermally  thick  interface  to  the  substrate. 


2 . 5 . 2 . 2  Boundary  Value  Solution 

Equation  (1)  was  solved  to  obtain  the  temperature  at  the  surface 
in  terms  of  the  flux  at  the  surface,  noting  that  the  substrate  is 
assumed  to  be  infinitely  thick,  and  so  the  solution  in  the 
substrate  contains  no  terms  propagating  toward  the  surface.  The 
solution  is: 


surface 

where 


A  C  C  +■  A  C  S  +  A  S,C„  +  A  S,S„ 
_cc_l_2 _ cs_  1  2  sc  12  ss  i  2 

B  C  C  +  B  C  S  „  +  B  S,C-B  S  ,  S  „ 
cc  1  2  cs  1  2  sc  12  ss!2 


(2) 


9 

CC 

‘6. 

m 

C9 

a2 

162 

m 

8C 

a 

3 

♦ 

‘8, 

m 

88 

a 

4 

+ 

164 

and 
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and 


B 

B 

B 

B 


cc 

cs 

sc 


ss 

where 


A 

A 

A 

A 


sc 

s  s 

cc 

cs 


1  +  h23X2 

*  g23h12x1//n 

h23X2//’r 

+  g23h12Xl//7 

h12xi//7 

g23 

h12Xl// 7 

+  2h12Xl  h  2  3X  2/lr 

e23  g12 

0 


g12  h23X2/'/* 


In  these  equations  Xj  ,  and  h^  are  dimensionless  parameters 

defined  as  follows: 


X1  '  > 


1/2 


x2  = 


1/2 


where 


c  .  I 


C  .  2 


.1 


f  .  are  characteristic  frequencies  above  which  thermal  diffusion 
c .  J 

wave  Interference  effects  are  absent,  due  to  the  critically 
damped  nature  of  oscillatory  solutions  to  the  diffusion  equation. 
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The  Modulation  frequency  f  Is  the  independent  variable  in  the 

problem . 


g12  * 

V5i 

g23 

gjj  are  ratios  of  effusivities 

of  the 

respective 

layers 

(  Ref 

When  glj  - 

1 ,  then  the  regions 

i  and 

j  are  thermally 

matched 

there  is  no  reflection  of  the 

thermal 

diffusion  wave 

a  t 

the 

interface . 

Thus,  the  ratio  of 

effusivities  is  to 

therma 1 

diffusion 

waves  at  an  interface  what 

the  ratio  of 

ref  ract i ve 

indices  is 

to  light  waves  at  an  interface. 

h12  = 

R12/R2 ,c 

ST 

to 

CO 

II 

R23/R3,c 

where 

R„ 

2  ,  c 

*  1/e2/lffc.l 

R 

3  ,  c 

=  1/e3,n'{c.2 

where 

6J  " 

y^j 

Here,  e.  is  the  thermal  effusivity  of  region  j,  and  R  .  is 
J  J  .  c 

characteristic  thermal  resistance  of  region  j,  for  a  layer  in 
region  J  whose  thickness  is  defined  by  the  characteristic 
frequency  of  layer  j  -  I  .  (Her.  (»),  h^  j  is  a  dimensionless 
parameter  proportional  to  the  contact  resistance  between  regions 
1  and  j,  but  also  with  imp  Mr  it  dependences  indicated  by  the 
definitions . 

In  equation  (2),  the  t  l  me  -  i  ndependen t  complex  phase  of  the 
heat  flux  at  the  surface  is  zero,  and  the  t ime - i ndependen t 
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complex  phase  of  the  temperature  at  the  surface  is  the  phase 
shift  between  the  heating  and  the  temperature  response.  This 
phase  Is  to  be  equated  to  the  phase  measured  in  the  experiment  by 
the  lock-in  analyzer  as  a  function  of  modulation  frequency  f. 


2. 5. 2. 3  Complex  Phase  of  Solution 

If  the  complex  components  of  Eq .  (2)  are  separated,  then  the 

complex  phase  is  given  by: 


i  -  f 


numerator 


denon ina  tor 


-  TT/4 


(3) 


where  numerator  and  denominator  refer  to  the  fraction  within  the 


parentheses,  and  the  v/4  is  due  to  the  divisor  C,.  In  this 
equation , 


6  .  =»  tan 

”  numerator 


f 


denominator 


=  tan 


-1 


I m ( numerator  )  1 


Re(numerator) 


Im( denominator ) 


R  e ( d  e nominator) 


where 

Re ( numerator )  =  a  ^(cCjCCg-sSjSSg 

♦  a  2 ( cCj  sc2 -ssj  cs2 

♦  a3(SClCC2'CSlSS2 

♦  a^sCjSCg-cSjCSg 

Im(  numerator  )  *  a  ^  (  ss ^  cc2  <-cc  t  ss2 
+  ClgtSSjSCg+CCjCSg 

+  a3(cslc<Vscl992 

♦  a4(cs1sc2+sc1cs2 


-  ( SS jCCg+CC jSSg) 

-  B2 ( ssj sc2 +cCj cs2  ) 
"  B3(csiCc2+sClss2) 

-  B4 (cs1 sc2 ^sCj cs2  ) 

*  BjIcc^Cg-sSjSSg) 

4  B2(CC19C2_SS1C92) 

♦  63(9C1CC2'C91992) 
+  B4 { SCj SCg-CSjCSg ) 
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and  where  the  corresponding  real  and  imaginary  parts  of  the 
denominator  are  obtained  by  letting 

1  ->  3 

2  ->  4 

3  ->  1 

4  ->  2 

in  the  and  coefficient  subscripts.  In  these  expressions, 
cCj  =  cosh(Xj)  cos(Xj) 
cSj  -  cosh(Xj)  sin(Xj) 
sCj  =  sinh(Xj)  cos(Xj) 
sSj  =  sinh(Xj)  sin(Xj) 

2.5.3  Graphical  Presentation  of  Theoretical  Solution 
The  theoretical  expression  for  phase  Eq.  (3)  contains  six 
adjustable  parameters  (g12>  g23 .  h12,  h23>  f£  ^  fc  2>  and  one 
independent  variable  f.  The  full  complexity  of  the  solution 
behavior  has  not  been  explored  yet.  One  illustrative  example  is 
presented  here.  Figure  5  shows  interference  oscillations  in  the 
phase  for  a  high  contrast  case  consisting  of  two  coatings  of  high 
thermal  conductivity  deposited  on  a  substrate  of  high  thermal 
conductivity,  and  with  Interfaces  of  low  thermal  conductivity 
between  the  coatings.  The  alternation  of  high  and  low 
conductivities  is  responsible  for  the  high  contrast  of  the 
interference  oscillations.  Points  to  be  noted  include  the 
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Phase  +  45  [degrees] 


Fig  5  Two-Layer  Coating  (High  Contrast  Case) 


following.  The  oscillations  damp  out  at  high  modulation 
frequencies  due  to  the  critically  damped  nature  of  thermal 
diffusion  waves.  The  frequencies  and  phases  at  the  four 
principal  phase  extrema  constitute  theoretical  data  structure  to 
be  matched  to  measured  data.  Note  that  these  features  span  six 
decades  of  modulation  frequencies.  The  present  measurement 
apparatus  will  tune  over  only  four  decades  (10  Hz  to  100  kHz). 

If  the  rear  surface  probe  can  be  applied  to  two-layer  coatings 
without  the  introduction  of  systematic  errors  (which  remains  to 
be  seen),  then  the  equivalent  of  at  least  one  independent  feature 
can  be  obtained  without  the  need  for  additional  dynamic  range  in 
modulation  frequency.  This  example  again  illustrates  the 
necessity  of  developing  the  greatest  possible  system  dynamic 
range  in  order  to  capture  data  over  a  sufficient  interval  of 
interference  oscillations  to  assure  uniqueness  in  fitting  the 
adjustable  parameters. 


3.0  CONCLUSION 

An  improved  measurement  system  for  performing  thermal 
transport  studies  In  coatings,  interfaces,  and  surfaces  has  been 
developed  and  implemented  ,\  major  source  of  systematic  error  in 
earlier  work  has  been  discovered  and  eliminated.  Measured  data 
are  in  qualitative  agreement  with  theoretically  predicted 
behavior;  numerical  fitting  of  the  data  constitutes  the  next 
research  phase.  Efforts  continue  to  characterize  in  detail 
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thermal  transport  in  simple,  thin  layered  structures.  Initial 
theoretical  modelling  of  more  complicated  two-layer  structures  is 
reported. 
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